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If an infinite flat elastic plate containing a circular holé is sub- 
jected to a loading which amounts to simple uniaxial tension at distances 
remote from the hole, the classical solution by Kirsch provides an evalua- 
tion of the stress components throughout the plate, — the loading 
| is sufficiently small. However if the tensile loading is progressively 
increased, there comes a point when Kirsch's solution becomes invalid, 
either through inelastic action at the most highly stressed regions in the 
plate, or through buckling of the plate from its original plane. The 
question of buckling under these circumstances has previously been dis- 
cussed only by Danis, who dealt experimentally with finite plates, and by 
Pellett who performed a theoretical study of an infinite plate. he 
present thesis pinpoints and corrects some errors in Pellett's ae 
and leads to the result that buckling impends when the tensile stress 
reaches the value 


S.- 77208 (t/a)? 


where E denotes Young's modulus and t/a denotes the ratio of || thick- 
ness to hole radius. This evaluation is for Poisson's ratio 7 '=0.3, 


a commonly used value, but evaluations are also made for other values of 


4 , indicating that the variation with respect to 2’ is quite small. 
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NOMENCLATURE 
a Hole radius 


A A set of unknown coefficients in a series describing 
plate deflection 


A An eigenvector 
b The half-width of an infinite strip 
B A coefficient matrix (Appendix B) 
Cc A coefficient matrix (ipoeudie C) 
d Hole diameter 
D Flexural rigidity of the plate = Eh?/12 (1 day?) 
E Young's Modulus 
h Ratio of plate thickness to hole radius, (t/a) 
T Identity matrix 
I A dummy integer index for the index n 
a I' A single subscript representing the double subscript IJ 
I, Integral in numerator of equation 2 
I, Integral in denominator of equation 2 
J A dummy integer index for the index m 
ax A single subscript representing the double subscript KL 
J Number of unknown coefficients = MN 
K A dummy integer index for the index n 
K Buckling coefficient 
L A dummy integer index for the index m 
é m An integer index 
M Number of terms in summation on m 
: ML Bending moment per unit length 
Mee Twisting moment per unit length 








NOMENCLATURE (Continued) 
An integer index 


Upper limit of summation on n . 


Membrane forces per unit length = a h, etc. 


A term in the quadruple series that expresses the 
bending energy integral 


The matrix of elements PLL 


A term in the quadruple series that expresses the 
membrane energy integral 


The matrix of elements Qs omer 


Distance from the hole center, also a coordinate 
of the cylindrical system 


A matrix of the eigenvectors of P 


Intermediate functions of N, M, andz (Appendix A) 
Externally applied stress 


The value ofgS for incipient buckling 


Cusemedlete Soambhina of N, M, andz (Appendix A) 
Plate thickness 

Component of displacement in the x direction 
Component of displacement in the y direction 


Component of displacement in the z direction (lateral deflection) 


Derivatives of plate deflection with respect to the variables 
indicated by subscripts 


A coordinate in the cartesian system 


A coefficient which is used to satisfy the boundary conditions 
for each term in n or m 





NOMENCLATURE (Continued) 


Xx An eigenvector (Appendix B) 
- y A coordinate of the cartesian system; also a scalar 
parameter in an eigenvalue problem (Appendix B) 
« XY A coefficient which is used to satisfy the boundary 
7 conditions for each term inn orm 
XY Vector of coefficients in Appendix C 
& The variation operator, (also see equation 20) 
€ The sum of the squares of errors (Appendix C) 
6 An angular coordinate in the cylindrical system 
us A scalar parameter in the eigenvalue problem 
v Poisson's ratio 


O45 % Fr Sy Normal stress components 


Gy, Ce Shear stress components 
ay A stress function (Appendix D) 
. OTHER NOTATION 


A double underline denotes a square matrix. (A) 


A single underline denotes a column vector. (A) 


T (as a superscript) denotes transposition of a matrix. (a*) 


-l (as a superscript) denotes inversion of a nonsingular matrix. at 
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Ne Introduction. 

The well known solution by Kirsch (6) for the membrane stresses in 
an infinite plate containing a circular hole and under a state of uni- 
axial tensile stress at an infinite distance shows that there are regions 
in the plate where one or both of the principal membrane stresses are 
compressive. Thus it seems reasonable to expect that for sufficiently 
thin plates under sufficiently great tensile stress, the plate may buckle 
out of its original plane. Danis (2) and Pellett (7) seem to be the only 
ones who have studied this problem. Danis performed tests on plates of 
finite dimensions and inferred, incorrectly we now believe, a certain type 
of ianicombanag of buckling stress upon ratio of plate width to hole diameter. 
Pellett considered the problem of the infinite plate, using an energy 
method to calculate the critical buckling stress. However, his analysis 
contains some errors which invalidate his results. These errors came to 
attention during attempts to deal with "Howland's problem", described in 
Appendix D hereof, and attention was inmediately focussed upon correcting 
Pellett's analysis and computing corrected results. Gratefyl acknowledge- 
ment is made to Pellett for his cooperation in verifying the nature of his 
errors and in offering suggestions for correcting the analysis. 

During the recomputations which were necessary, using the corrected 
theory, two additional improvements were made on what Pellett had done 
earlier. For one thing, the nature of the analysis requires that the 
buckling stress be inferred as a lower bound of a number of determinations 
each of which is itself slightly in error. In the present analysis a 
least-squares inference is made so as to maximize the objectivity of the 
inference, Also, although it might be expected that variation in Poisson's 
ratio would have but slight influence upon the buckling stress, a quanti- 


tative comparison is of interest, and this is accomplished herein. 
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Figure 1 


The specific problem under investigation is illustrated in Figure 
ol The center'of the circular hole is taken as the origin of a set of 
cylindrical coordinates (r, @, z) and a set of cartesian coordinates (x, 
y, 2). The z-axis, not shown, forms a right-handed orthogonal set with 
the x and y axes. The radius of the hole is @ (See Figure 1) and the = 
plate thickness is t. For convenience we temporarily take the unit of 
length to be one hole radius, so that the corresponding hole radius is 
unity and the plate thickness is h = t/a. However, the influence of 
variations in hole radius is properly re-introduced in reporting final 
results. 

At a very large distance from the hole it is assumed that the state 
of stress is a simple uniaxial tensile stress S in the x direction, If 
the value of s is sufficiently small, Kirsch's solution correctly des- 
cribes the stress ddineetbucion in the entire plate. However, if S is 
iaeeaused , Kirech"e solution will eventually become invalid, either through 
inelastic action of engineering materials stressed above their proportion- 


al limits, or through buckling of the plate out of its original plane. 
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Kirsch's solution, along with a knowledge of the properties and behavior 
of the material, is sufficient to permit determining the value of S for 
which inelastic action begins. The problem of interest herein is deter- 
mining the value of S for which buckling is imminent, under the presump- 
tion that the material is still perfectly elastic. 

This problem will be called the problem of incipient buckling. The 
critical value of S is that above which Kirsch's solution is no longer 
valid because of buckling. However, no attempt ig made herein to dis- 
cuss the post-buckling problem, that is, the problem of determining stress 
and deflection distributions for values of S greater than that required 
for incipient buckling. It is clear, however, that information about this 
evidently difficult problem would be of value. 

Similarly, since all plates in existence are finite in all dimensions, 
it would be of interest to consider finite plates. The infinite plate is 
an evident good case with which to start because, besides offering fewer 
mathematical difficulties than the infinitely long strip of finite width 
(hereafter called the infinite strip) or the finite plate, there is avail- 
able a tractable solution for membrane stresses. Before the errors in 
Pellett's analysis came to attention,the writer was attempting to deal 
with the problem of the infinite strip. Although no useful results were 
obtained for this much more difficult problem, nevertheless some details 
of the attempts are reported in Appendix.D for whatever use they may be to 
others who will work on this problem. 

It is proper to delineate the relation of this thesis to that of 
Pellett. The viewpoints are essentially the same. The fundamental equa- 
tions from which the analysis proceeds are essentially the same except that 
Pellett made an error in his derivation which amounts to a factor of two. 


In the present thesis no derivation is made, the appropriate equations being 
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taken directly from a standard reference (8). Most of the remainder of 
the analysis is essentially the same in the two theses, except for differ- 
ences in notation, a slight simplification of equations occurring in Ap- 


pendix A, and the early introduction, in the present thesis, of matrix 


that the present thesis does not include certain tabular (and similar) 
listings atid presentations given by Pellett, the reason being that if a 

; ati ' 
reader requires these details and is unable to reconstruct them easily 
for himself, he may find them in Pellett's thesis. Furthermore, as 
mentioned earlier, the present thesis contains an evaluation of the influ- 
ence of variations in Poisson's ratio and a least-squares extrapolation 
for the lowest beiticnl stress (see Appendix C), neither of which appears 


in Pellett's work. 


notation. Many of the details in the developments are so nearly identical a 
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2. Analytical Procedure. 
Timoshenko and Gere (8) discuss two methods of solving problems such 
: as the one posed here. One method is to seek solutions of the equilibrium 


equation * 


tere B+ Sy = SN de an, B+2N ze | (1) 

y % OK Y oy RY Ou dy 

| where NaN > and iy are membrane forces, w is the lateral deflection, 

| and D is the plate flexural rigidity, D = Eh/12(1- 7 eh. This equation is 
to be solved for w, subject to certain boundary conditions on w. 

Timoshenko and Gere also discuss what they gall the energy method. 

Because of the evident mathematical difficulties which would be encountered 
in attempting to use the first method for dealing with the present buckling 
problem, it was decided to use the second method, that is, the energy method 


which is contained in the following equation: (equation 9-3 of reference 


a 8) ** 
w E (= 
- Sone ef 4 pao ) J aay 
where X is a common ercoiee \ N,N, and N_, so that 
x sy xy 
ei Ny can ee i= 10, 


* 
Both equations 1 and 2 assume that the membrane forces N._, N , and 
are independent of the buckled deflection. This assumption is’ satisfied 
oy the problem of incipient buckling. 


eve 
Pellett (7) derived this expression independently; however, he made 
an error which, together with another error (see page 24 ), invalidated 
. his results. 
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are the membrane forces per unit length. 


For the calculation of ¥ it is necessary to find an expression for 


1 


w which satisfies the boundary conditions and minimizes ¥ , i.e., the 


variation of ¥ (in equation 2) must be zero. 


aw 


Ww 
dopting the io = = : - . 
Adopting t notation We se My so? etc., and trans 


forming this equation into cylindrical coordinates, with the proper limits 


of integration, results in the equation 


Lyrik d, = ay 


oe soften OF oT r 
+ a-aN[ viel aes Me) = Ct - es] arae 
and i Pw we +N, wa +2N. “| rdrde (5) 
‘ Gal a 


The values of Nt N'y and Lar which are to be used in equation 5 


are obtained from Kirsch's solution (6) for the membrane stress distri- 


bution, namely, 


CO, = ${c- 4) + (+3, - 4) cos 26| 


(6) 
e+ $Lliy A) ~ (14 Br) eo 20] m 


- ae 

The equation (eq 5) which expresses Tg in this paper is the negative 
of that which Timoshenko and Gere use. This accounts for equation (3) 
being slightly different than that given by Timoshenko and Gere. 
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where S is the applied tensile stress (see Figure 1). 
The membrane forces may be expressed as NL = oh Ny = Oph. and 
= Cag It is now apparent that ¥ may be interpreted as S$ in the 


6 


above equations. Using the notation orn = eae » etc., equations 3 


and 5 become respectively 


I, +9 1,+=0 (9) 


aT «wo 


2 hf ft 6iwi+ Gq + +27. ‘ewe | rede de (10) 


The assumed expression for w must have sufficient disposable para- 
meters that the boundary conditions may be satisfied and that the actual 
lateral deflection may be closely approximated, Pellett (7) chose the 


following excellent expression for w 


oo wo a 
wine) ea > [1 + Yap gl {, =f \ May 


By) aie) 4 sin me 


where Aum and Bm are arbitrary coefficients that are varied in order to 
minimize S in equation 9, coefficients Xm and Yum having been chosen so 
as to satisfy the boundary conditions. 

From Figure 1 it is seen that there is loading symmetry with respect 
to both the x and y axes. This means that when w is written as a Fourier 
series in 9 (eq. 11) the series contains either cosine or sine terms, but 
not both. Additionally, the requirement of symmetry around both the x and 
y axes means that only odd or even terms in m need be used. Thus the 


Fourier portion of equation 11 reduces to: 
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(a) Ann cos me, Mi O52, gestae s (Sym A)* 
Symmetry with respect to the x and y axes 
(b) or Aum Sosen@cully Me vlc. By -5 phe in (Sym B) 
Symmetry with respect to the x axis 
Antisymmetry with respect to the y axis 
(c) or Amn Sin m@> jm =id, 8, Sy ea. (Sym C) 
Symmetry with respect to the y axis 
Antisymmetry with respect to the x axis 
(d) or Aun sin nO, Mt S705 4, anes (Sym D) 
Antisymnetry with respect to the x and y axes 
This results in four separate expressions for w, one for each possible 
condition of symmetry. The full series for w could be used, but if this 
is done a large atari of a storage would be wasted because in any 
solution three quarters of the coefficients, 4am and Bm? would be zero. 
It will later be shown that each of these expressions for w represents a 
different buckling mode. 
The function w must also satisfy certain boundary conditions. The 
F . re 
stress or strain distribution is that resulting from the Kirsch solution 
(boundary conditions for which are already satisfied) and the infinitesimal 


lateral deflections w. 


The conditions which must be satisfied at r = -<O , namely 


Wiz We We =O (12) 


* 
These symmetry designations are used here simply to distinguish between 
the different types of symmetry. The letters agree with ascending order of 
buckling stress, Sor? obtained in the analysis which follows. 
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are obviously satisfied by the choice of w given by equation 11 since w 
and all its derivatives vanish for r = -o, 

The edge of the hole is free. Following Timoshenko and Woinowsky- 
Krieger (9) and using the cylindrical coordinate system established in 


Figure 1, the conditions for a free edge are 


ic, -¥ SE | = 0 (3) 
r=\ 
[mM], = 6 a 


While the work necessary to show that the expression for w, equation 11, 
satisfies these equations is not difficult, it is somewhat involved and 
therefore is shown separately in Appendix A. Farmulas for the quantities 
Qw Mee and M.. are also given in Appendix A. 

It is sufficient for our purposes here to say that unique values of 
Xm and Yom are found which allow each term of equation 11 to satisfy the 
free edge boundary conditions for the values of n and m chosen. 

The procedure is now straightforward even if very complicated and 
detailed. Inasmuch as Pellett's thesis containg the details (particularly 
his Appendices B and C which cannot be abbreviated without losing intel- 
ligibility) this material will not be repeated here. The reader who wishes 
to examine the details more closely is referred to Pellett's thesis (7). 
Briefly the procedure is as follows. The Kirsch solution, equations 6 — 
8, with the stress § Eastowal out, is substituted into equation 10. The 
appropriate partial derivatives of w are formed from equation 11 and sub- 
stituted into equations 4 and 10. Note that in these integrals the de- 
rivatives of w appear to the second degree. Since w itself and each of 
its derivatives is represented as a double series, such second degree 
terms are quadruple series. Evidently there is considerable manipulation 


to be done and a great deal of "book-keeping" which must be managed. 
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The actual integrations are elementary although Pellett finds it conven- 
ient to devote an entire appendix (his Appendix C) to some of the details. 
The upshot is that after the substitutions are made, the multiplications 


performed, and the integrations carried out, the energy integrals I, and 


I, take the form . 
190% 4 % 4, Prot Ars AR <uease 
To=bh 3 3% 4 Oe Arg Age te 

where PLL and ORL are numerical coefficients. Actually, as may be 

seen from the ating LaMocne for I, and I,» they are really functions of 


Poisson's ratio, 3/ , but because of evident difficulties no attempt has 
been made to retain the literal dependence upon 7 ; rather, the evalua- 
tions which are actually performed incorporate a numerical value for ¥ 
at an early stage. Most of the calculations were actually carried out with zy 
a value of 7 = 0.3. ‘Rover, supplementary calculations were performed 
using other values of 7/ so as to assess the dependence of the final re- 
sults upon the choice of) 7 . 

In sguneiows 15 and 16, the quantities Ary and Avy, represent the 
double series of as-yet unknown coefficients in equation 11. It does not 
seem at all possible to carry out the remainder of the evaluation in literal 
terms so that recourse is made to a digital computer evaluation. Obviously, 
chen, it is necessary to truncate the series given in equation 11. Thus 


we take the summations in equation 11 so as to include N terms, M terms 


*The unknown coefficients for the cosine series, Anm, only will be 
used henceforth, although these could be the sine series coefficients, , 
Bam equally as well. 
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respectively. Later it will be seen that different choices of N and M 
lead to slightly different values for the criti¢al stress, and advantage 
é will be taken of this fact to extrapolate the results so as to obtain a 
better evaluation than would otherwise be passible within the limits of 
computer storage and time availability. 
With such a truncation, it is clear that the number of unknown 
quantities A is J = N-M, and it is convenient to order these quantities 


into a single, singly subscripted array, which we do as follows. We write 


ie TJ 
where 
T'=(I1-1)M+JI+1 
Using the same device, the quadruply subscripted quantities PLKL and 
i! QF KL are written as doubly subscripted arrays Piyt and Qrrze where 


each of the subscripts runs from 1 to J. 

To make this operation clearer, the following example is offered. 
Suppose that the value M = 4 has been chasen (the value of N does not 
matter in this renumbering) and we wish to convert 85432 = 2(5,4) (3,2) to 


a doubly subscripted quantity. We have I = 5, J = 4, K = 3, and L = 2. 


From the formula above, we calculate I' = 21, J' 11. Thus we write 
the quadruply subscripted quantity 25439 (i.e., que sub five four three 
two) as Q514,11 (i.e., que sub twenty-one eleven), With this change in 
notation, equations 15 and 16 may be written as 


I= D 2 Five! om aN gg (17) 


rat sel 





Le Wot > A oA. 8) 


vey T's 


A further simplification results if matrix notation is adopted at 
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this time. In matrix form equations 17 and 18 become 


1 ORAL, .. tare el 


where A ete LA, ‘", As} 


and P, Q are symmetrical matrices of order (IxJ). The unknown A;'s must 


be such that S is stationary. (a minimum) for arbitrary variations in the 


SA = col { 54,,+- +, SAS} (20) 
a = 1. ._ @ @ eal (21) 
“OA OAy yg » As 


we proceed as follows* 


$5902 ee SA 
iL, 


_— 


£ (+s BT 


re) 


(22) 


Using equation 19 and the analysis in reference (5) 


A SA 2 


* 
The’ second error in Pellett's work occurred when he differentiated 
incorrectly in the process of obtaining the variation of S. 
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Since §&S is to'be zero for arbitrary §A, we are led to 


2| a al deo (24) 


2 
When the substitution “\= ss = eaiae is made, the following 


matrix eigenvalue problem results 


[>Q+P]A=0 (25) 


The method used to solve this matrix eigenvalue problem is shown in 
Appendix D. There are as many solutions, eigenvalues », » as there are 
elements of A. The smallest positive eigenvalue, which corresponds to 
the minimum applied tensile stress, is of the most interest. Larger 
positive values of “. correspond to more complex buckled mode shapes 
while the negative values correspond to buckled modes due to applied 
compressive stress. 

The value of the buckling stress is recovered from the definition 
of a . That is, we have 

a 

5. = fafeli-w)]e* 
However, we recognize, from consideration of the dimensions involved, that 
we can restore the analysis to a case of hole radius other than unity 
simply by writing t/a in place of h. Also, we denote the term in square 
brackets by the symbol K, which we call the buckling coefficient. Note 
that the effect of varying Poisson's ratio cannot be inferred from the 
expression in square brackets inasmuch as Poisson's ratio has been inti- 


mately involved in the entire analysis, starting with equations 1 and 2. 


Ze 





To investigate the effect of varying Poisson's ratio it is necessary to re- 
5 j PFs. 1 


calculate the quantities P and Qa KL using different values of 2’ * 


IJKL 
and also to vary 7/ in the expression (1 - wv * 


Thus we finally arrive at the expression 
Ss = KE (t/a)? (26) 
cr 


where K is a function only of z/ . Different modes, that is, different 
values of the eigenvalues ‘\ , and different assumed conditions of 
symmetry, lead to different evaluations for K, of which, naturally, the 


smallest is of greatest interest by far. 


*Different values of z/ can be introduced in the computations by 
changing card 0011 in the main program and card 0003 in the subroutine 
BOUND; see computer program listing, Appendix E. 


26 





3. Results. 

Computed values of K in equation 26 are shown in Table 2; these 
values assume 7/ = 0.3. N and M in Table 2 are the number of terms 
taken in n and m in the assumed deflection, equation 11. As both N and 
M are increased, the value of K appears to approach a limit. In Figures 
2 and 3, K is plotted versus M for varying values of N. In Figure 4, K 
is plotted versus N for varying values of M. The trend toward a limit 
is more obviously shown in these curves. 

Because of limitations on computer storage space and time, N = 8 
and M = 10 are the practical upper limitsfor the number of terms that 
can be used in the series expression for deflection. 

Similar curves and tables could be exhibited for other values of 7. 
However, since they are similar to those for z/ = 0.3, they are omitted. 


A tabulation of K for N = 8 and M = 10 is shown in Table 1 below. 


Table l 
Values of the Buckling Coefficient, K, for N = 8 
and M = 10 
Mode Sym.* w= 0.0 v= 0.2 y= 0.3 w= 0.4 v= 0.5 
ii A 1.786 1.723 1.726 1.752 1.808 
2 B 2.106 2.026 2.025 2.051 2.109 
3 Cc 3.332 
4 D 3.651 
5 A 5.646 
6 B 6.377 
7 Cc 7.851 


* 
4 These symmetry types correspond to the symmetries discussed on page 
0 e " 
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Table 2 
Values of the Buckling Coefficient, K, for Various Values of M and N (7 = 0.3) 


Saw ~ ee et ae 


1.7450 = 1.7269 1.7267 
- L.#352 1.7254 
= 1.7348 


Mode 1 








2S SS. 2a Se a a 








3.3349 31.3331 3.3325 3.3323 
SeaSZ1 


: 2.0255 2.0254 
: 2.0252 2.0250 
es) 







3.3348 Jadsae 3.3324 






be Se Be ee ee a 


3.6587 = 2 
- 3.6529 3.6524 3.6522 
- 3.6518 3.6513 3.6512 


(no values were calculated where blanks are indicated) 
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Since the lowest value of K and corresponding critical stress is of 
primary interest and 7 = 0.3 is a common value, more time was devoted 
analyzing the value of K for wz = 0.3, mode l. 

Figures 2 and 4 indicate that the value of K approaches a lower limit 
or minimum as the number of terms in the expression for deflection (equa- 
tion 11) is increased. When the deflection is approximated by truncating os 
the series in equation 11 we are introducing hidden unknown constraints 
on the deflection. These constraints tend to stiffen the plate and re- 
quire a greater applied load to cause buckling. As we more nearly ap- 
proach the true deflection, by increasing the number of terms used in the 
series expression for deflection, we more closely approach the value of K 
for the actual plate. It is in this sense, of approaching the true value 
of K, that minimum is used. 

While it appears from Figures 2 and 4 that K is approaching a limit- 
ing value, it is not apparent exactly what that value is. It appears that 
the minimum value of K is about 1.72. Accordingly, a method was sought 
by which these curves could be extrapolated very objectively in order to 
obtain a better minimum value of K. 

Using a least squares method* of curve fitting the minimum value of 
K was determined to be 1.720. 

For any symmetry type such as Sym A (symmetric with respect to both 


the x and y axes) there are several** values of K which correspond to 


* 
See Appendix C for the details involved in this technique. » 


ek 
An infinite number of such values would be expected if-.it had been 
possible to use an infinite number of unknown coefficients 3} however 


A 
with the use of truncated series there are only a finite number of such 
values. 








increasingly complex buckled deflections, each of which retains the same 
symmetry. For example, for Sym A., with 7/ = 0.3, N= 8, and M = 10, 
the following values of K were determined, 

Mode = 1 5 9 12 


K 1.726 5.646 9.205 12.086 


Figures 5 and 6 illustrate the buckled deflection for modes 1 and 2 
respectively. The more complex deflection of mode 5 is illustrated in 
Figure 7. It can also be seen that the symmetry is the same as for mode 
1, 

Figure 8 is taken from the experimental work of Danis (2). His 
results for finite plates are plotted as K versus d/b, the ratio of hole 
diameter to plate width. The buckling coefficient for mode 1 is plotted 
on the y axis (d/b = 0). 

There appears to be some agreement between the value of buckling co- 
efficient calculated in this thesis and Danis' experimental work, The 
minimum value of K calculated was 1.720 (for 7 = 0.3), while the lower 
values that Danis found were between 1.5 and 2.5. For small d/b ratios 
Danis found values of K that were much higher than this. A possible 
explanation may be that the strain gages he used to measure the buckling 
reinforced the areas near the edge of the hole so as to resist buckling. 
If the area covered by the strain gage was to become greater in propor- 
tion to the area that was likely to buckle, the reinforcing effect would 
become greater. This is exactly the situation for Danis’ low values of 
d/b, as he used smaller holes in the same size plates to vary the ratio 
of d/b, This reinforcing effect would raise the critical stress required 


to cause buckling. Also, it seems reasonable to assume that the effect 


of the edge of the plate would not be much different for d/b = 0.3 or 
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. Figure 8 - Comparison with Experimental Data of Danis 
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d/b = 0.1. In Figure 8 the solid curves are Danis'‘s and the dashed curve 
is a postulation of what the results should be, on the basis of the above 
discussion. 

As seen in Table 1, Poisson's ratio has a small effect on the buckling 
coefficient. This is also illustrated in Figure 9, where K is plotted - 
versus Poisson's ratio for mode 1 with N = 8 andM = 10. The difference 
between K for 2/7 = 0.3 and 7/= 0.0 is less than 4%. 

This small effect of y/ may be helpful in the solution of other 
buckling problems. Consider some other buckling problem which has not 
been successfully solved by any method of attack, but for which there is 
reason to believe that the effect of z/ will be small. If the problem 
is simplified by letting 2/ = 0, it may become amenable to attack, per- 
haps by a theoretical analysis. There would be good reason to conjec- 
ture that the solution for 2 = 0 would be an acceptable solution for - 


practical values of z7/ such as 7 = Oe 
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4. Conclusions. 

The critical applied tensile stress that will cause incipient 
buckling in an infinite plate pierced by a circular hole can be calcu- 
lated using energy methods. The minimum value’ of critical stress, Peel 
occurs when the buckled deflection is aiuimetric about the x and y axes 


and is equal to 
2 
S 21.720 E(t/a) 
cr 


where E is Young's modulus of elasticity, t is the plate thickness, and 
ais the hole radius. This evaluation is for Poisson's ratio zw = 0.3, 
a commonly used value. 

The effect of Poisson's ratio on the critical stress is not 
great. The maximum difference observed for euiuee of z/ from 0 to 


@.5 is less than 5%. 
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APPENDIX A* 
Satisfaction of the Boundary Conditions 
The free edge boundary conditions for the hole are stated in 


equations 13 and 14 as 


(Q, -+ Spt) =O (A-1) 
(M,) > O07 (A-2) 


fal 


Timoshenko and Woinowsky-Krieger (9) show the expressions for Q- 


M_., and MM. to be 


rt 
2. 2 
Q,=-D&[ ea +te +h Syl 3) 
> 2 wi ie D om 
M, = -D[ Pe +o( poe +4 Be) | (A-4) 
: mos JA 
Mi = DU-wLt Be - BT] (4-5) 


P 
Adopting the notation = = we etc., and performing the appropriate 


differentiations and substitutions, the free edge boundary conditions be- 








come 
2-2 B-y) _ 
| Were ST ‘er a “ 2 r2 Wree { =¥ Woo] = . (A-6) 
W, _ 
[ Wr + (m+ e)| =O (A-7) 


* 

Appendices A and B are essentially the work of Pellett (1). Only 
minor changes in notation have been made to correspond to changes in the 
main text. 
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The lateral plate deflection w, as expressed in equation 11 is 


me ao * 33 nontge m@ 
winds > > fit Se ce 4 E (A-8) 


Bim sin mO 


where the oprtisl ante Ain and B am are varied to minimize the plate 
energy, and the coefticiente: Xm and oe are chosen so that the boundary 
conditions expressed by equatloes A-6 and A-7 are satisfied. 

The derivatives of w are tabulated in Table 3. When these deriva- 
tives are substituted into equations A-6 and A-7, and evaluated at r = 1, 


the following equations result 


am Cos mo 
> S$ 5,4) Kin wien | = 0 (A-9) 


n=p M30 Bam sin mO 


= Anm cos me 
S (+ tka, wor oe = 0 (A-10) 


Bam sin m@ 


5, =(nei)(nti¢h-v) - my (A~11) 


; a ‘ 
4, =(3-p)m +(2-v)mi(nsi) —(n+i) (n+i+2) (A-12) 
for i=0, 1, 2 
Equations A-9 and A-10 must hold for all values of 9, which can be 


true only if 


+ = - 
85 . 4 — “2 Teen ° (A-13) 
£ +0 +e Y = 0 A- 
fe) 1 mm 2 nm = 


for any value of n and m. Thus, a value of Xm and Ym can be found for 


each value of n and m, such that the boundary conditions are satisfied. 
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The results of solving equations A-13 and A-14 simultaneously are 


Xm (82 ty ~ 89 to) / (8, ty - 8, ty) (A-15) 


¥ = (s5 t, - sy ty) / (s) ty - 85 t,) (A~16) 
The above results are incorporated into SUBROUTINE BOUND (see Ap- 


pendix E) which generates values of Xam and Yam for given values of n 


and m. 
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w' = pS = [rvs a ae eas I vie nm cog me = 


nH=1 m20 
Bim sin 


aot) aie eer nm cos mO 
We ==> > [ne +(nti)K ‘y (n42) > aad {™ 

oo a(n) aw cos mO 
Wer - Pad { n(n-i) y «+ (nti\(n4 2) Yu. an * + (n42) (n4 3) a 4 (os Bey , 


nat ry 
Ps Bren run mo 


> S sie nA “ Kam cos mO 
vers ES foro e™ samara AMrameaneny ee] (mer 
Pa net om=0 
am sth m& 


=D eo e 
. nahn F m +3 am sin m 
Wee - 7 > [a r ) + (nan) ee v mt + (n42) meta) Ws ‘ 
n=l oO 
Ws qos ‘me 
= = —(n +2) ~(n+3 nm cos me 
Wree ~ [n _ + nal he F mw (nt 2) Dt ee ee | m sO: \ 


ta sin me 


-ina') | nn aes mO 
i ows YY * Yous SS : 
Bam sin m& 


-n -tn+1) =(n¥2) m= Nam sin me 
Y + X.Y + Yam ™ 
cog mO 


Table 3 
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8 

Vw 
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7M 3 
aie 3M g 
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Be 
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Tabulation of the Derivatives of w 





APPENDIX B*¥ 
Solution of the Eigenvalue Problem (Q \ + a & = 2 


The canonical form of the matrix eigenvalue problem is 
x - Bg = 9 (B-1) 


where I is the identity matrix, X an unknown eigenvector, x is an un- 
known scalar parameter (eigenvalue), and B a known square coefficient 
matrix. 


Equation B-1 has a nontrivial solution if and only if the determinant 


det (ZL x - B) = 0 (B-2) 
If this determinant is expanded, a polynomial equation in x results. How- 
ever, if the determinant is of order greater than three or four, it is 
generally inefficient to obtain this polynomial equation and then solve it. 
In the actual numerical calculation in this thesis,the so-called Jacobi 
method (3) was employed. 

Consider the eigenvalue problem (equation 25) 

QA +P A = 0 (B-3) 
where both Q and P are symmetric matrices. Dividing equation B-3 by » 
and making the substitution 

y= - z (B-4) 
gives 


®y -Q@ A= 0 (B-5) 


Any symmetric matrix with real elements can be reduced to a product 


of the eigenvector and eigenvalue matrices, so that 


* 
See footnote page 42. - 
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e- hee ) (B-6) 
where D is a diagonal matrix of the eigenvalues (d,) of P, and R is a 
matrix of normalized eigenvectors (x,) obtained from a solution of the 
problem 


(L 4, -P = 0 (B-7) 


i DE 
Note that R= x since the matrix P is symmetric and has orthogonal 
eigenvectors. 


' Equation B-5 may now be written 


(RDR'y - Q) A = 0 (B-8) 


Premultiplying by R’, we get 


(Ry - BQ) A = 0 (B-9) 
Factoring R out of A = IA = RRA, and letting 

Y= RA (B-10) 
the problem reduces to 

wy, - R™QR) Y = 0 (B-11) 


Since P is a result of the bending energy integral it can be 
expected to be positive definite*. D, the matrix of eigenvalues of P, 
therefore contains all positive terms and the square roots of its elements 


are real numbers; (only the positive square roots will be used). Pre- 


multiplying by p-, factoring, and making the substitution 
xe at 
Z=) (B-12) 
equation B-11 reduces to . 
(ly - p~#x"gpp”*) Z = 0 (B-13) 


* 
During the numerical computations it was found that P is in fact 
positive definite. 
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and becomes 
(ly = B) 2 = 2 (B-14) 
which is of canonical form, where 


% (B-15) 


B= D*R QRD 
% -% / , 

Note that D“ and D “ are easily obtained by taking the square roots, . 
and the square root reciprocals, respectively, of the diagonal terms. 
Also note that B is symmetric. The problem is now well suited for solu- 
tion by use of the Jacobi method. It ig seen that inaccuracies in obtain- 
ing eigenvectors R would cause compounded errors in the final solution 
for A and . 


The problem is returned to original form with the substitutions 


aN 


il 
' 
<I 


(B-16) 
and 


“5 (B-17) ‘ 


A=R D4 
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APPENDIX C 
Least Squares Extrapolation of K 
Since it is desirable to estimate the minimum buckling coefficient, 
the following method of extrapolation was used. 
It is assumed that the buckling coefficient K is a function of the 


number of terms used in the series for w. Such a function may be 


Y Y Y Y Y 
Key t+ 3+ e+ t+eeesy, en 
M M N N nM 


where K is a calculated buckling coefficient, N and M are the number of 
terms in the series for deflection, and Y,> Yooreee¥e are unknown coef- 
ficients. 

If the coefficients of equation C-1 can be determined then ¥y will 
be an estimate of the value of K as both N and M approach @ ., 

For a number of calculated values of K and known values of N and M 


a series of equations of the above form could be written in matrix form 


as follows 


bs to ae (c-2) 
where C2 Ci Cio ity K= K, 
Cor &o0 Yo Ry 
Cc, %6 Ry 


Z is the number of values of K that were calculated 


Cia =1, C = 1M,”, etc. 


A common method of fitting an equation to data points is the so- 
called method of least squares. In this method the sum of the squares of 


the errors, resulting from the use of a proposed equation such as equation 
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C-1, is minimized. That is, the variation of the sum of the squares of the 
errors with respect to the unknown coefficients, in this case the Yi» is 
set equal to zero. 

Using © to denote the sum of the squares of the errors, the result- 


ing equation is 


e=(¥-*) (XY - DW (c-3) 
This equation may be expanded to the form 
do age t - tee ek & OE (C~4) 


If Y has variation § XY, € will have variation* 


&& = = gy = 2(Y c¢- KC) bY (c+5) 
and for this to vanish for arbitrary variation in Y, it is sufficient 
that 
v= Kg gy” (6-6) : 
It is easy to see that this gives a minimum since 
we ar 
sy = 2G (c-7) 


A digital computer was used to evaluate equation C-6 for the 15 
values of K shown in Table 1, for mode 1, and zw = 0.3. The computer 


program used is shown in Appendix E. The results obtained are shown be- 


low 
yy = 1.720 
Yo = 0.1479 
Y3 = 30.93 
Y, = 0.1202 
Y, = 3-553 F 
Y¢6 = 12.17 ; 





4 
Differentiation of expressions of this type is discussed by 
Frazer, Duncan and Collar (4), 
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Other functions expressing the dependence of K upon N and M could 
be assumed, such as using odd negative exponents for both N and M. 
When this was done the coefficients associated with the odd negative 
exponents were found to be either zero or much smaller in magnitude 


than the coefficients listed above. 
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APPENDIX D 
Howland's Problem 

It was the author's original purpose, before discovering the errors 
in Pellett's analysis, to extend this work to the case of the circular 
hole in an infinite strip having finite width and having no loadings on the 
finite boundaries. The first requirement for a solution to this problem is 
to determine the membrane stress distribution. 

Since Howland (5) has discussed this problem and provided an itera. 
tive solution for the membrane stress distribution, it is appropriate to 
call this problem by his name, Howland's problem, However, his solution 
was not considered very tractable for use in an energy approach to the 
buckling problem. Howland's solution involves the sum of a finite number 
of terms, d aa each of which is an Airy stress function, to approximate the 
Airy stress function for the problem. In theory each of the co) -e contains 
an infinite number of terms; however, in practice only a finite number of 
terms can be used. Howland's solution becomes invalid at relatively small 
distances from the hole (approximately one plate width or less). This dif- 
ficulty could be eliminated by using a larger number of O's and more terms 
in each of the > ‘s- However, his method of determining the coefficients is 
quite involved. For these reasons, several attempts were made at solving 
Howland's problem by other methods. 

Neither of the two approaches to this problem which are presented in 
this appendix was successful. It is hoped that the information presented 
here may prove useful to someone else who may be interested in the solution 


of this problem, 
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Description of the Problem 


Consider a strip of finite width, 2b, infinitely extended in the 
° other direction, pierced by a hole at the center, and subjected to uniform 


tensile stress, S, at infinity as shown in figure (10). 


y 
- 
ate S 


Figure 10 





It is convenient to employ two sets of coordinates: x and y are cartesian 
coordinates and r and 9 are cylindrical coordinates. The boundaries in 
the y direction are at +b. Without loss of generality all dimensions can 
be considered in units of hole radius. This reduces all quantities to 

. dimensionless ratios with the hole radius equal to unity. 


The boundary conditions for this problem are as follows: 


(o,),., =0 (D-1) 
¢ cma =o (D-2) 
( Oy \yetp =O (D-3) 
( Tghe gh =O (D-4) 
6, jie 3 (D-5) 
a (D-6) 


In each of the attacks on the problem which will be described below, 
. the point of departure is assuming a stress function from which the stress 
components may be determined,as with an Airy stress function, so as to 
satisfy the laws of equilibrium. However, instead of requiring this 


function to satisfy the biharmonic equation, it is selected so as to 
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boundary conditions are satisfied by including corresponding terms, multi- 
plied by undetermined factors (Lagrange multipliers) in the strain energy e 
function. Budiansky and Hu (1) present a detailed example of this method 
in minimizing the strain energy in a certain plate problem and at the same 
time satisfying the boundary conditions. 
In the first approach to this problem, consider a stress function 
such that the boundary conditions on the hole and the edges of the strip 
are satisfied. The Lagrangian multiplier is to be used to satisfy the 


boundary conditions at x = c©® . Such a function might be as follows: 


O(r,6) =f (r,e)+ 4lr,e)* h(r,e) (D-7) 
where f (r,@) is such that the boundary conditions at the edge of the plate 
are satisfied, and g (r,@) is such that the boundary conditions on the 
hole are satisfied. The function h(r,@) is to be general enough that the 
actual stress distribution in the plate can be closely approximated. The 


functions selected were 


(1 - Esinte =(1-¥) (D-8) 


f 


i 


file 4Y (0-9) 


N . 
2 >. A om r cos2mO (D-10) 


nzo0 Meo 


SS ee 
uf {t 


It can be shown that the boundary conditions D-1, 2, 3, 4, are in 
fact satisfied by the above choice of the functions f and g. The function 


minimize the membrane strain energy in the plate, a number of disposable 

parameters having been introduced so as to permit minimization and so as to 

permit a stress distribution approaching the correct distribution. The 

h appears to be quite general, considering that positive exponents of r 





must be excluded in order for the resultant stresses to remain finite as x 


(and r) approach infinity. 
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The boundary conditions at x = sO , namely equations D-5 and D-6, 
are yet to be satisfied. In the limit as x approaches infinity oc. ap- 


proaches C.. and r approaches x. The boundary condition D-5 becomes 





(Gh aes = 5 (D-11) 
This results in a relatively simple equation in terms of the stress 


funetion > A 


ae 
4 > Age =5 (D-12) 
m=O 
The difficulty in the solution arises from the fact that the stress 
at x = c© is dependent upon y. The stress varies in a parabolic manner 


as shown in Figure ll. 





Figure 11 


When the stress is integrated from y = -b to y = +b in order to 
determine the average value of Cw the result is zero, i.e., there is no 
net loading on the strip, just local stresses at x = + oe. 

It is interesting to note that while the above expression for d 
appears quite complicated, that it can be reduced to a very compact form. 
If the functions f and g are expressed as binomial expansions and cos 2m@ 
is expressed as powers of sin 9, equations D-8, D-9, and D-10 become 


respective 


- 


{ 
(TY sre) 0-13) 





9 = Bal) (-r ia 


n2O m26 Ks 


N tm 
h = » a het Cae Es sin (D-15) : 
° 


2 = m2 Cael L/axdiim=%)l 


wig, SESS S 


1=O jw n=O mea K=O 


and Des) Pile - (?) (2\ eh 


The resulting expression for > is 


(5) + >» , 2 K+2i 
> P Z Des) — ar: one? 


The form of equation D-16 considerably simplifies the job of dif- 
ferentiating in order to determine the stresses. It also simplifies 
the expression for the strain energy to a great extent. This form is also 
well suited for use in a computer solution to the problem. 

Howland's method of solution involves starting with a stress func- 
tion d, » which results in the Kirsch solution for stress. This would 
result in Sy and buy stresses on the boundaries at y = +b. He then adds 
a stress function fo) such that the stress on the edges is reduced to 
zero but there would have to be loadings on the boundary of the hole. 
Howland then adds another stress function 5: which reduces the stresses 
at the hole to zero but introduces loadings at the edges of the plate. 
However, the stresses at the edge due to >,+ ?, + @ are different 
than those due to b S alone. He continues this process until the stresses 
on the edges in the neighborhood of the hole are reduced to. an acceptably 


low value. 
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Howland's solution is of the form 


P= &+h+o,+  --- +, pay 


where each fa) ons of the form 


b. = 28, ee3 = ( B. r+ Ge aa Cos 2nO (D-18) 


nel 


and each Pe. is of the form 


n=O 


Rather than use Howland's method of solving for the coefficients of each 


of the > "s the following stress function is proposed. 


| $= h +644 0-0 


where > _ is the same as Howland's p " b 1 and d, are of the follow- 


ing form 


of 


> (D. ene EL at, cos 2ne (D-21) 


hzO 


b 
®, 


This > is equivalent to Howland's vo) 4 


tt 


oo , ~on+2 
Sg (Br = C, r ik ) cosane tAlar (D-22) 


he) 


The intention is to solve for 
the unknown coefficients using energy methods rather than Howland's alter- 


x Nating procedure. Since d _ Satisfies the boundary condition at x equal 


to infinity, the boundary conditions of the sum d 5 + d 9 are 


( ox hues *~. (ae =O 


In order to Satisfy these conditions all of the coefficients Dd. and 
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EA must be equal to zero, i.e., > 1 =.0. 
Since d a satisfies the conditions for zero tractiona on the edge 


of the hole, p 2 must also satisfy these conditions. In order for 


2 2 
to satisfy these conditions the coefficient A and all the coefficients 
Bo and Cc. must be equal to zero, i.e., o, = 0. = 
Hence, b = o °° which we know does not satisfy the boundary condi- 


tions at y = +b. 


Suggestions for solving Howland's Problem 


Since as a practical matter it is necessary to truncate the series 
used in Howland's solution, his solution provides valid results only for 
a limited range of values of x. Accordingly, the so-called "finite element 
method" applied to a finite plate should provide as good a solution to 
Howland's problem as Howland himself was able to provide. 
One is motivated to deal with the problem of the plate of finite 
width because all plates in actual practice have finite width. Inasmuch 3 
as all real plates also have finite length, it can be convincingly argued 
that the finite element method applied to a finite plate would produce 
more meaningful results than would a solution to Howland's problem, how- 
ever the latter might be achieved. 
Those familiar with finite element techniques contend that the finite 
element method can be employed to solve the buckling problem directly. 
Also, with what would appear to be developments within the state of the 
art, it would be possible not only to deal with the incipient buckling 
problem, but also to deal with the post-buckling problem in which the ‘ 


membrane stress distribution is coupled with the lateral deflection. 
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APPENDIX E 
Computer Program Listings 

The computer program listings on the following pages were used in 
the solution of the problem. All programs are written in FORTRAN IV for 
use on an IBM System 360 computer. 

The program was originally written by Pellett (7); however, modifica- 
tions were made to the program to reflect changes in the solution to the 
problem as noted previously. 

The plate thickness, N, and M are read into the computer by data card 
(see comment cards regarding input requirements). Young's modulus is 
established by card 0010 in the main program. Poisson's ratio is estab- 
lished by card 0011 in the main program and card 0003 in subroutine BOUND. 

The output is in the form of a column vector EIVU which contains 
the values of Sur in psi for the parameters of plate thickness, Poisson's 
ratio, Young's modulus, M, and N. This-is followed by a deflection 


contour map for the lowest value of tensile buckling stress. 
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A 2 
a. = 1.720 E (t/a) 


where E denotes Young's modulus and t/a denotes the ratio of plate thickness to 
hole radius. This evaluation is for Poisson's ratio wv = 0.3, a commonly 
used value, but evaluations are also made for other values of z’ , indicating 
that the variation with respect to z’ is quite small. 
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